Analytical relation between entropy production and 
quantum Lyapunov exponents for Gaussian 
bipartite systems 

K. M. Fonseca Romero 2 ' 3 , Julia E. Parreira 1 , L. A. M. Souza 1 , 
M. C. Nemes 1 ' 3 , and W. Wreszinski 3 

1 Departamento de Fi'sica, Institute de Ciencias Exatas, Universidade Federal de 
Minas Gerais, CP 702, CEP 30161-970, Belo Horizonte, Minas Gerais, Brazil; 
2 Departamento de Fi'sica, Facultad de Ciencias, Universidad Nacional, Ciudad 
Universitaria, Bogota, Colombia; 

institute de Fi'sica, Universidade de Sao Paulo, CP 66318 CEP 05389-970, Sao 
Paulo, SP, Brazil. 

E-mail: carolina@fisica.ufmg.br 

Abstract. We study and compare the information loss of a large class of Gaussian 
bipartite systems. It includes the usual Caldeira-Leggett type model as well as Anosov 
models (parametric oscillators, the inverted oscillator environment, etc), which exhibit 
instability, one of the most important characteristics of chaotic systems. We establish a 
rigorous connection between the quantum Lyapunov exponents and coherence loss and 
show that in the case of unstable environments, coherence loss is completely determined 
by the upper quantum Lyapunov exponent, a behavior which is more universal than 
that of the Caldeira-Leggett type model. 



PACS numbers: 05.45.Mt,03.65.Ud,05.70.Ln,03.67.-a 



Submitted to: J. Phys. A: Math. Gen. 



Analytical relation between entropy production and quantum Lyapunov. 



2 



1. Introduction. 

Application of Quantum Mechanics to real world many-body systems meets with several 
difficulties, both of conceptual and pragmatic nature. A quantum mechanical system, 
which consists of at least two interacting subsystems exhibits a completely nonclassical 
property called "entanglement" [lj. It is a nonclassical correlation between systems 
which exists even between well separated subsystems [2j |3]. This unique property of 
a quantum system is nowadays viewed as a powerful resource in quantum information 
theory and quantum computation [4J. Also, on the conceptual side, "entanglement" 
with an environment (degrees of freedom which inevitably interact with the quantum 
system of interest) is the key mechanism to explain why typically quantum effects are 
not observed in macroscopic systems [5j [6]. Therefore the rate at which pure initial 
quantum states loose their potentiality to retain information represents at the same 
time an old and a very modern problem. 

The first theoretical implementation of the system plus environment dynamics, 
in both the weak and strong damping regimes, was proposed and carried out by 
Caldeira and Leggett [7]. They represent the environment as an infinity of oscillators, 
weakly coupled to the system of interest. Several successful descriptions of experiments, 
specifically in quantum optics, were based on this type of model. More recently questions 
have been raised about the role of classical chaos on the information loss process for 
coupled quantum systems. In particular, the role of chaos in the decoherence process 
[H] has become a matter of active research and also a matter of debate P, QUI [TT] . 

One of the important issues now is: which type of environment is the more 
effective: the one involving an infinite number of degrees of freedom (quantum brownian 
motion environment QBME), as proposed by Caldeira and Leggett, or an environment 
consisting of a single degree of freedom which presents one or more characteristics of 
chaotic behavior? 

In the standard model of decoherence, a joint system consisting of a system (S) 
coupled to an environment (E) is described by a Hamiltonian 

Hj = H s + H E + V, (1.1) 

where H$ is the system's Hamiltonian, He the environment's, and V the interaction 
Hamiltonian between S and E. Let, for simplicity, the state space Tis of the system be 
two-dimensional, and 

H s = Xa z , (1.2) 

with a z the Pauli matrix with eigenvalues +1 and —1 corresponding to spin "up" and 
"down" along the z-axis, and respective eigenfunctions ip + and ?/>_. In the coherent 
superposition 

ip = a+V>+ + a -if)-, (1.3) 

with probability amplitudes a + and a_ such that \a + \ 2 + |«_| 2 = 1, o~ z has no definite 
value. The environment state (Hilbert) space is TCg, and we suppose E is initially in the 
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state ^s(O) G 'He an d the joint system in the pure (nonentangled) state 

^j(O) = V ® V»b(o) g W 5 ® (1.4) 
We also assume for simplicity that the interaction in (1.1) is of the form 

V = ^zQe- (1.5) 

Then, after a time t, the joint system will be in the state (units are chosen such that 
Ti=l) 

ipjit) = (a + exp[— iXt]i/; + ®i^%{t) + a_ exp^'At]^- <8> '4 ) e{^))i (1.6a) 

with 

^(t) = exp[-i(H E ± hQeMe{0). (1.6b) 

We thus expect that the sensitivity to perturbation ±[iQe of a chaotic environment 
Hamiltonian He will, by (1.6a) and (1.6b), result in quick decoherence: 

\(ipUt),M-t))\ <cexp[-dt] (1.7) 

where d and c are positive constants. The exponential "orbital" instability of He is 
expected to lead via (1.6b) to the exponential decay of the "overlaps" (1.7) (think of 
Gaussians centered at the corresponding orbits): the decoherence rate should therefore 
be proportional to the rate at which the environment is able to explore its phase 
space. These considerations are meant to explain, intuitively, why a (large) collection 
of harmonic oscillators - being stable systems - could, in principle, be less effective for 
purposes of decoherence than one sole unstable system. 

Recently, R. Blume-Kohout and W. H. Zurek [12] analyzed decoherence due to a toy 
model for the environment, an inverted harmonic oscillator environment (IHOE) - see 
also [131 HH E3 where conclusions similar to ours were obtained in chaotic systems. This 
is a nontrivial analytical tractable model whose importance as realized in [12J is that 
it shares one very important characteristic with chaotic systems, namely exponential 
instability, defined by a positive upper Lyapunov exponent [IE]. Systems of this sort 
are called Anosov systems. We shall be rather dealing with generalized Anosov systems, 
in the sense that the flows take place in a manifold which is not compact. A simple 
example is the flow in the (p, q) plane defined (with A positive) by -£ = Xp and -S = —Xq. 
The q coordinates of any two points moving with the flow get closer and closer as time 
proceeds, but the p coordinates separate exponentially fast, and hence the two points 
move apart exponentially fast [IT] . The directions of stable and unstable manifolds may 
vary from point to point, however, in contrast to the above example: this is the case of 
the parametric oscillator. The quantum analogues of the classical Anosov systems have 
been dubbed quantum Anosov systems |18) : they exhibit a dynamic behavior which 
is at the same time rich and universal and they are characterized by having a positive 
upper quantum Lyapunov exponent, as defined in [19] and in the forthcoming (11751) . and 
will be the subject of the present paper. 
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Given the importance of the subject, the difficulty in obtaining mathematically 
sound exact results which are eventually able to shed light onto these complex questions, 
we feel that a thorough analytical investigation of the QBME vs. IHOE and its 
generalizations to more realistic models of the class of linearly coupled quantum Anosov 
systems, can be extremely useful. We focus on the limited but important class of 
gaussian states. In this context, a rigorous demonstration of the growth of the system's 
reduced von Neumann entropy with the upper quantum Lyapunov exponent for the 
aforementioned class of systems is given, showing thus that the rate of information 
loss is of more universal nature than that for the usual quantum brownian motion 
environment. Our analytical results are only possible due to the fact that for (possibly 
nonunitary) gaussian dynamics the rate of information loss is completely governed 
by a combination of quadratures, i.e. the covariance matrix, also called Schrodinger 
generalized uncertainty principle. 

2. Relation between information loss and the covariance matrix for 
gaussian states. 



We start by recalling a well known result [20J [2T] . The most general 1-D gaussian state 
can be written as 



p G = V(a)S( r ,4 > )p u S\r,4 > )V\a), (1) 



where 



T>(a) = exp (cra^ — a*a*j , (2) 

S(r,<P)= exp g (e^ 2 - e^a 2 )) (3) 

are the displacement and the squeezing operator respectively and p u the thermal density 
operator with average number of excitations z/, 

p u = exp fin ( ^ d^a\ . (4) 



1 + v V + 

Note that except in the improper limit of v — > 0, which corresponds to a squeezed 
state, density matrix ([!]) is not a pure state. Its Schroodinger determinant is given by 



D= ( °vv *«, \ ={y + \f (5) 
V a qp a m ) 2 

where a xy = | tr (p {x, y})—tr(xp) tr(yp) and x, y are either the position or its canonicaly 
conjugate momentum. We moreover have for the (von Neumann) entropy 

S[p) = -tr (plnp) = (z/ + 1) ln(z/ + 1) - z/lnz/. (6) 

Comparing ([6]) with (J5]) we obtain the relation between the entropy and the Schroodinger 
determinant 

S= i(v / ^+l)ln(v^D + l)-i(v / ^-l)ln(v / ^-l))-ln2. (7) 
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It is worthwhile to point out that holds for a wide range of (gaussian) states, 
ranging from thermal (with S(r,(p) = 1, T>(a) = 1 in ([T])) to generalized coherent states 
(the limit v — > in with S = 1 in (CQ)); ([7]) is thus a fundamental relation for 
gaussian states. 



3. Quantum Brownian Motion Environment (QBME). 

We next consider an example of QBME widely used in the context of Quantum Optics, 
whose dynamics (in the Born-Markov approximation) is governed by the Liouvillian 



£ = — iu 



a) a, 



a) — a) a • — • a) a 



+ k(n B + 1) (2a 

+ kfiB (2a) • a — aa) • — • ad^ . (8) 

The entropy is a function of u(t) only (see ([6])), where 

K*) = vVata(t)) 2 -^tat(tKa(i) ~\ (9) 

and 

<M.(*) = (KO) + ~) cosh(2r(0))e- 2fct + (n B + ~)(1 - e - 2fci ), (10) 

<W*yaa(t) = e_4fct ((^(°) + \) 2 sinh 2 (2r(0))) . (11) 

Notice that the entropy for this type of system may be obtained without resorting to 
the calculation of the full density matrix, due to the fact that the information exchange 
between the system and bath is completely determined by a simple combination of 
quadratures. Also the entropy saturates for long enough times. Moreover, from equation 
(J9]) it is apparent that information loss depends both on the environment constants as 
well as on the initial conditions. 



4. Inverted Harmonic Environment (IHE). 

We now revisit the inverted harmonic oscillator example. Following (1.1), we have now 

H s = pj/2 + cufxl/2 (12) 
H E = pl/2 + uj 2 2 xl/2 (13) 
V =Axix 2 (14) 

The above Hamiltonian models two coupled harmonic oscillators: the first one, is a 
regular oscillator, with unit frequency, describing the system; the second one, an inverted 
oscillator with squared frequency uj\ = —A 2 , describing the environment. This is a 
nontrivial tractable model. Its importance, as has been realized in reference p2], relies 
on its relation to chaos. Chaotic behavior requires besides stretching and folding in phase 
space, instabilities characterized by positive Lyapunov exponents. Although the above 
hamiltonian is defined on the whole of phase space and does not exhibit chaotic behavior, 
it displays positive Lyapunov exponents like all Anosov systems. In this section we 
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establish a connection between Lyapunov exponents and the rate of information loss, in 
this model, which we generalize to periodic quadratic Anosov systems in the next section 
and exemplify with a physically more sound system of coupled parametric oscillators. 
The solution of the Heisenberg equations of motion for the vector z = x 2 ,P2) T is 

as follows [22] 

z(t) = Ge^G^z, (15) 

where G and L are constant matrices. The former is invertible and the latter is diagonal 
with elements ±A = ±A(A,a) and ±fi = ±iB(A,a), where A(A, a) and B(A,a) are 
real and positive constants. The IHE Hamiltonian is a degenerate case of Floquet's 
theorem, see the following. Notice the existence of an unstable mode, with positive 
classical Lyapunov exponent A. In reference [19] . upper quantum Lyapunov exponents 
were defined, for systems of one degree of freedom with momentum p and position x, as 
follows 

A = sup limsup -In \\[L a , A(t, to)]\\ (16) 

aeR 2 t-KX t 

were a = (a p ,a x ), and A(t,t ) is any bounded operator, i.e. with finite norm given by 
||y4|| = sup^, || J A'0||/||'0||, where if) belongs to the Hilbert space L 2 (R, dx), and evolved 
in the Heisenberg picture. This norm is nothing but the natural extension of the usual 
matrix norm on the square integrable sequences to an infinite-dimensional Hilbert space. 
The operator L a = a p p + a x x is the generator of phase-space translations along a 
direction a, and thus ( 1T61) has a close resemblance to the definition of the classical 
Lyapunov exponents (see, e.g., [23]), but exploits the unitary nature of the quantum 
dynamics (see [19] for a discussion of these points). 

Also in reference [19] the definition (fT6l) was applied for one of the simplest 
paradigms of the transition from regular to unstable behavior in classical mechanics, 
namely the parametrically driven oscillator, given by 

H(t) = \p 2 + \f{t)x\ (17) 

where f(t) is a periodic function. The dynamical solution for x(t, 0) and p(t, 0) is a 
formally one dimensional version of ffTBT) . but with G being a periodic matrix with 
the same period of the driving and L a traceless matrix: this is Floquet's theorem 
and the eigenvalues of L are known as Floquet exponents (see [19] for the proof and 
references). Applying definition (TIB"]) it was found that the upper quantum Lyapunov 
exponent is precisely the real part of the maximal Floquet exponent, which is positive 
in the classical instability region of the parametric oscillator, and zero in the stability 
region. The generalization to two degrees of freedom is straightforward [23] and, applied 
to the IHE shows that the upper quantum Lyapunov exponent (fTBl) is A (see equation 
eingenvalues of the IHE). Let now S r = — tr (pi lnpi), the von Neumann entropy of the 
reduced density matrix, with pi = tr 2 p, with p of the form ([I])-®, but now with v — > 0, 
i.e., the initial state being a tensor product of generalized coherent states. Since pi is 
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also Gaussian, it continues to satisfy d7j). In both cases, IHE and parametric oscillator 
(j!7p . the determinant is found to be of the form 

D(t)= £ C ab e^ x+b ^\ (18) 

a,b=-2 

where the constant coefficients C a b are given by cumbersome not particularly 
enlightening expressions. These coefficients, which are zero unless the sum of a and b is 
an even number, depend on the variances of the initial state and will be omitted here. 
For instance, one can find C a b by looking at the asymptotic behavior of (EJ). For large 
enough times t we have D(t) ~ C 2 oe 2A * and the reduced entropy can be approximated 
by S r ~ ln(-/D - 1/2) + 1 ~ lnD/2. This behavior, 

S r (t) ~ ln(C 20 )/2 + Xt (19) 

is precisely the one conjectured by Zurek several years ago [12]. Since p\ is also gaussian, 
(J7J) applies to p\ and thus the linear growth of S r (t) is interpreted, in the present case 
of a bipartite system, as a coherence loss, which, as we see, depends only on the 
upper quantum Lyapunov exponent (uqLe) and is thus universal for this class of 
systems, in contrast to QBME. Note that complexity in the present models (they are 
not explicitly soluble classically and are a paradigm of the transition from regular to 
irregular behavior, see [25]) is brought about by the external field, which depends on 
time in a nonlinear way! For the models with frequency varying almost periodically with 
time - included in the present treatment - see [26]: it may be considered as a prototype 
of complexity. This is one of the reasons why our extension of [12] is significant: Zurek's 
IHOE is not complex! 

As remarked in [12], section II, p. 032104-2, the IHE Hamiltonian is unphysical 
in that in general the directions of stable and unstable manifolds can vary from point 
to point. Thus it is of special interest to be able to consider as system 2 an arbitrary 
quantum Anosov system. In what follows we illustrate this by taking both systems of 
the form (fTTl) . a case in which most computations can be done explicitly. This model 
is particularly rich: there is even a transition from a stability to an instability region, 
which has a physical interpretation in a model of quadrupole radio-frequency traps 
(Paul- Penning traps) (see [19] and references given there). 

For the sake of comparison, in figure 1 we plot the entropy obtained for a quantum 
Brownian motion environment (BME) and for a IHO reservoir (parameters: (i) for the 
BME u = 1, r = 1, k = 0.5, n B = 10, u = 1, (j) = 0; (ii) for the IHO ln(C 20 )/2 = 1.4 
and A = 1). Note that the increasing in the entropy in the case of BME saturates when 
the system thermalizes with the reservoir, while the increase in the case of IHO is linear 
(in the asymptotic regime). 
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Figure 1. Time evolution of the entropy for a quantum Brownian motion environment 
(solid curve) and for a inverted Harmonic Oscillator reservoir (dash curve). 



5. Coupled parametric oscillators. 

Consider now a more realistic example, physically relevant in connection with Paul traps 
|27J, which models two coupled parametric oscillators with hamiltonian 

H(t) = H s (t) + H E (t) + V (20) 



where 



and 



Hs{t) = -pl + (w? - gcos2t) xj, 



H E{t) = - f 2 + (wf - gcos2tj x\ 



V = |(xi - x 2 ) 2 . 



(21) 
(22) 

(23) 



In what follows we assume g > and lo\ > uj\. We remark that fl20]) corresponds 
to the standard model (1.1). The solution of the Heisenberg equations of motion for 
the position and momentum operators can be written in terms of Mathieu functions. 
For example, the momentum of the second oscillator, p2{t) is given by the expression 
ui(t)xi(0) + U2(t)pi(0) + U3(t)x2(0) + W4(t)j3 2 (0) , where the functions u a (t), a = 1, 2, 3, 4 
are given, respectively, by 



Mi = sin(2#) 



S(a-,q,0) . v S{a+,q,0) . . 

-C(a_, q, t) — C(a+, q, t) 



u 2 = cos 



m 3 = sin(26>) 



2D 2 

, S(a_,g,0) 



D 2 

C(a-,q,0) 
2Do 



2Di 

C (a_, g, t) — sin 6* 



C(a+,q,t) 



S(a-, q, t) — S(a+, q, t) 



2D X 



2 C(a-,q,0) . , . 2 C(a+,g,0) 
m 4 = cos 6* — i(a_,g,t)-sm t) 



£> 2 



S(a + ,q,t). 



(24) 

(25) 
(26) 
(27) 
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The parameters a± = Ul 2 ^ 2 + g ± y g 2 + 2 ^ depend on the constants appearing in 
the Hamiltonian. The angle 9 appearing in the expressions above is determined from 
the equality: tan(2#) = , 2 . The functions C(a,q,t) and S(a,q,t) are the usual 

Mathieu cosine and sine functions, and C(a,q,t) and S(a,q,t) their time derivatives. 
From the expression C(a,q,t) = \ (e' l ^ t P(a,q,t) + e~^ t P(a,q,—t) S j , where P(a,q,t) 
is a periodic function with period tt and = <p(a, q) is the so-called characteristic 
exponent, it can be seen that, if 3(0) 7^ the solutions exhibit an unstable behavior. 
Otherwise the solutions are uniformly bounded. We are interested here in a situation 
in which the second uncoupled forced oscillator (system) is in the stability region and 
the other (environment) is in the unstable region. From the solutions of the Heisenberg 
equations of motion and the definition of the upper quantum Lyapunov exponent A, it is 
a simple matter to show that A = \^(<f)(ai, q))\. The connection with the information loss 
is made through the evaluation of Schroodinger determinant for the second oscillator, 
which turns out to be of the form of equation ( 1T81) . In this case, however, the coefficients 
C a b are periodic functions of time. Moreover A = |S(0(ai,g))| and /i = i(p(a 2 ,q). For 
an initial gaussian state, S(t n ) ~ ln(C*2o)/2 + Xt n , where t n = 2nn. The coefficients C ao 
are obtained from the asymptotic behavior of the entropy ([6]). 

6. Conclusions. 

In spite of exhibiting a complex, rich dynamical behavior, bipartite open quantum 
Anosov systems display, in their instability region of parameters, a reduced von 
Neumann entropy with linear growth determined by their uqLe, a universal behavior 
for this class of systems, in contrast to QBME. This also permitted us to find a unified 
description of two opposite regimes, viz. of small and large coupling. It may, however, 
be argued that (|T8|) is induced by classical mechanics: the value of the uqLe for the 
present models equals the classical maximal Floquet exponent, which is positive in 
the region where the classical parameter E — ^ f f(t) dt lies in the instability region 
("gap") of Hill's equation + f{t)x = [28J, as may abstracted from the last section. 
This is a consequence of the fact that the dynamics of the observables is classical. The 
space of states is, nevertheless, subject to the laws of quantum mechanics, without 
classical analogue, even in these cases: for this reason, for gaussian initial states, 
the entropy growth f|T9l is to be interpreted, due to ([7]), as coherence loss, a purely 
quantum mechanical concept. Indeed, in contrast to quantum mechanical entropy, the 
classical entropy is not an adequate quantity in information theory, because it is not even 
positive in general: ([29], Prop.l). Negative entropies are due to the fact that classical 
density distributions may be concentrated in regions of phase space < h, contradicting 
the uncertainty principle ( [30] , [31]). Modified "semiclassical entropies" without this 
inconvenience have been defined by Wehrl [31], but they incorporate the Hilbert space 
structure of quantum physics, i.e, of a countable discrete basis of states. Entanglement 
has indeed been described in terms of this modified entropy [32]. Thus an explanation 
by "classical entanglement" of phenomena of information loss is not possible in a full 
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classical treatment involving the entropy. 

A still different approach is to regard entanglement in terms of the Wigner function, 
i.e., in phase space, which illuminates quite different and surprising aspects [33J. 

We should caution the reader that, even in the framework of (time-dependent) 
quadratic Hamiltonians, our approach misses an important element which is present in 
fully chaotic systems: folding, which is characteristic of a compact phase space. Indeed, 
other papers have reported on the growth of relative entropy in the context of quantized 
chaotic systems ([31], [35]). In contrast to the present treatment they rely, however, 
on approximations (based on random matrix theory in [31], semiclassical analysis in 
[35J), numerical studies and certain special constructions (the nonunitary step in [35] 
is constructed to mimic the presence of diffusion) and special concepts (a particularly 
phase-space measure of complexity was used in [34]). We hope that the discrete Weyl- 
Wigner formalism developed in [36] may be used to extend our treatment to some 
quantum chaotic systems. 

As a final remark, we mention that other parameters could be invoked for 
environments characterized by generic chaotic systems: the diffusion coefficient is one of 
them, which is conjectured to be related to the Lyapunov exponent and the Hausdorff 
dimension |37j . 
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